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Abstract 

We study the case of brane world models with an additional Gauss-Bonnet term 
in the presence of a bulk scalar field which interacts non-minimally with gravity, via 
a possible interaction term of the form —^R<p 2 . The Einstein equations and the 
junction conditions on the brane are formulated, in the case of the bulk scalar field. 
Static solutions of this model are obtained by solving numerically the Einstein equations 
with the appropriate boundary conditions on the brane. Finally, we present graphically 
and comment these solutions for several values of the free parameters of the model. 

1 Introduction 

Hidden extra dimensions is a way to extend four dimensional conventional physics. Accord- 
ing to the Kaluza Klein picture, extra dimensions are compactified to a very small length 
scale (naturally the Planck scale), and as a result space time appears to be effectively four 
dimensional, insofar as low energies are concerned. An alternative picture, which explains 
why extra dimensions are hidden, is the brane world scenario. In this case all ordinary 
particles are assumed to be localized on a three-dimensional brane (our world) which is em- 
bedded in a multidimensional manifold (bulk), while gravitons are free to propagate in the 
bulk [U [2] . The brane world scenario has attracted great attention, as it predicts interesting 
phenomenology even at the TeV scale [3], and puts on a new basis fundamental problems 
such as the hierarchy and the cosmological constant problem. According to the nature of 
the extra dimensions, there are two standard brane world models with gravity: a) in the 
case of flat extra dimensions we have the ADD-model [3], and b) in the case of warped extra 
dimensions we have the first and second Randall-Shundrum models jl]. 

Brane Models with bulk scalar fields have also been considered. In Ref. j5] the existence 
of a bulk scalar field serves as a way to resolve the cosmological constant problem (see also j6] 
and references therein) . However, the solutions of the model suffer from naked singularities, 
and this implies a fine-tunning for the five dimensional cosmological constant. An alternative 
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role of the scalar field is presented in Ref. [7]. In this case the scalar field creates the 
brane dynamically, by forming a kink topological defect toward the extra dimension. Also, 
nontopological defects appear for a special periodic form of the potential [8] . 

In Refs. [9j [10] we have studied brane world models with a bulk scalar field non-minimally 
coupled with gravity, via a possible interaction term of the form —l/2^R<f) 2 . Our motivation 
was to construct brane models which incorporate a layer-phase mechanism [TTl [T2"l [T51 [T4 J for 
the localization of the ordinary particles on the brane. Similar models for phenomenological 
reasons have also been consider in [T5| [T6] . 

Particularly, in Ref. [9J, by solving numerically the Einstein equations in the case of a 
non-minimally coupled scalar field, we have obtained three classes of static solutions with 
different features, in appropriate regions of the free parameters of the model. The numerical 
approach we apply in [9J is suitable for the general case of the potential. However, we have 
focused to the case of a massless scalar field assuming the simplest generally accepted form 
of the potential: V((f>) = A0 4 . 

An analytical study of the same brane model, with a nonminimally coupled scalar, can be 
found in Ref. JT7J. In this paper analytical solutions have been found by choosing suitably 
the potential for the scalar field. For cosmological implications of these models with the 
nonminimal coupling you can see |18j . 

Furthermore, 5D models in the framework of Brans-Dicke theory have been also studied. 
Particularly, in Ref. [15] . static analytical solutions are constructed for a special class of 
potentials, and the stabilization of the extra dimension is discussed. 

An alternative way to go beyond the Randall-Sundrum model, is to add a Gauss-Bonnet 
term to the gravity action, which, in the case of D > 5, has the property to keep equations of 
motion to second order. In the case of D = 4 the Gauss-Bonnet term is topological and does 
not contribute to the classical equations of motion, however it has nontrivial consequences at 
the level of the conserved currents of the theory [20] • It is worth noting that, a Gauss-Bonnet 
term arises as the leading order quantum correction in the case of the heterotic string theory. 
The Randall-Sundrum model with an additional Gauss-Bonnet term has been considered in 
[2T| 1221 [23l I21] . Brane models with scalar fields and Gauss-Bonnet gravity has been also 
examined in [251 123 [271 28J , while for cosmological implications see for example [29l 130] . 

In this paper we will study brane world models with Gauss-Bonnet gravity in the presence 
of a non minimally coupled scalar field. The motivation is to examine how Gauss-Bonnet 
gravity affects the three classes of solutions we obtained in our previous work [9j. 

2 Gauss-Bonnet gravity, brane world models, and non- 
minimal coupling 

In this section we will study brane world models with an action of the form: 




(1) 



where d 5 x = d 4 xdz, and z parameterizes the extra dimension. 
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As we see the above mentioned lagrangian consists of three terms. The first term corre- 
sponds to the lagrangian of the RS2-model: 



I ^1 \ni hrane )\ 

C<s = \/\u\ | 2^ " A " r- I < 2 ! 

where A is the five-dimensional cosmological constant and a is the brane tension. In addition 
R is the five- dimensional Ricci scalar, g is the determinant of the five- dimensional metric 
tensor g^N {M,N = 0, 1, ...,4), and g brane is the determinant of the induced metric on the 
brane. We adopt the mostly plus sign convention [3T] . 

The second term is the well-known lagrangian for Gauss-Bonnet gravity: 

r — ( fi> 2 A £? R M7V _L R r>MNKL\ fo\ 

l-GB — - [ft — ^tiMNti + ti MNKL ti J (6) 

where we have introduced a new coupling b in conection with the Gauss-Bonnet term. 
Assuming that this term arises from string theory (tree level) the coupling b takes only 
positive values. However, in this parer we will consider that the range of b covers all the 
real axis. 

The scalar field part of the lagrangian is 

^ = ^\[-\g MN Vu<l>VN^-\i<l> 2 R-V^)) (4) 

where the potential is assumed to be of the standard form V(<f>) = A0 4 . 
The Einstein equations, which correspond to the action of Eq. (1) are 



b J\ g (brane)\ 

G M n - t; h mn + A 9mn + a5(z)- — = — g^M^N = t mn ( 5 ) 
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where the energy momentum tensor for the scalar field is 
Tmn = Va/^V^ - (?a./v[^ PE Vp0V s + V{4>)\ - 2^V M V 7 v0 + 2^g MN U<f> + ^ 2 G MN (6) 
and Hmn is the Lanczos tensor 

Hmn = C-gb 9mn — ^RRmn + <$RmlR L n + &RmknlR KL — ^RmklpRn KLP (?) 
The equation of motion for the scalar field is 

n<f>-^<f>R-V'((j>) = Q (8) 

The above equation is not independent from the Einstein equations (5), as it is equivalent 
to the conservation equation V M T$ N = 0, where Tj$ N is given by Eq. (6). 
We are looking for static solutions of the form 

ds 2 = a 2 (z)(—dxl + dx\ + dx\ + dx\) + dz 2 , <\> = (f>(z) (9) 
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From the Einstein Equations (Eq. (5)) we obtain two independent equations: 

Gu - h -H u +Aga + a5(z)g u = I# } (10) 
G zz -^H zz + Ag zz =T^ (11) 

where i = 0,1, 2, 3. 

If we set a{z) = e A ^ and use Eqs. (9), (10) and (11) we get 

F 1 = : 3(1 - i<p 2 {z)) (A"(z) + 2A'(z) 2 ) - 12bA'(z)* - \2bA \zf A" \z) + A (12) 
+ (i - 2^\zf + V{cj>{z)) - 2^(z)</>"(z) - ^A\z)<P(z)<P\z) + a8{z) = 

F 2 = : 6(1 - ^\z))A\z) 2 - \2bA\zf + A - ~0'(z) 2 + V{<f>) - 8^ / ( 2 )0(^)0 , (z) = (13) 

From Eq. (9) for the scalar field we obtain 

F 3 = : -<f>"{z) - AA\z)<P\z) - £ (bA"(z) + 20A'(z) 2 ) <j>(z) + V\<P) = (14) 

Note that Eq. (14) can be found from the combination — AA'(z)(F 1 — F 2 ) + F 2 — of 
Eqs. (12) (F 1 = 0) and (13) (F 2 = 0). As Eqs. (12), (13) and (14) are not independent, the 
static solutions of the model can be obtained by integrating numerically the second order 
differential equations (12) and (14). The first order differential equation (13) is an integral 
of motion of Eq. (12) and (14), and acts as a constraint between A(z), (p(z) and their first 
derivatives. 

In particular we choose to solve numerically the following second order differential equa- 
tions: 

3(1 - ^ 2 (z))A"(z) - \2bA'{z) 2 A"{z) + (1 - 2^'(z) 2 - 2^{z)<P"{z) 
+ 2£A'(z)(t>(z)4>'(z) + a8(z) = (15) 

- (j)"(z) - 4A'(z)(f)'(z) - f (8A"(z) + 20A'(z) 2 ) <f>(z) + V'{<f>) = (16) 

Eq. (15) is obtained as the difference between Eqs. (12) and (13). 

This system of differential equations is quite complicated, so we will not look for analytical 
solutions. On the other hand we will try to solve it numerically. For the numerical integration 
(z > 0) of Eqs. (15) and (16) it is necessary to know the values of A(0), 0(0), A'(0 + ) and 
</>'(0 + ). These values are determined by the junction conditions (see Eqs. (17) and (18) 
below) and the constraint Eq. (13). 

The junction conditions on the brane are obtained from Eq. (15) and (16) 

6(1 - £0 2 (O))A'(O + ) - 8bA'{0 + ) 3 - 4^(0)^(0+) + a = (17) 
0'(O + ) + 8£A'(O + )0(O) = (18) 

if we take into account that A'(zfA"(z) = l/3(A'(z) 3 )'. For details see pH [251 M, [22] . 

The values of 0(0), A'(0 + ) and 0'(O + ) can be found by solving numerically the system 
of algebraic equations (17), (18) and equation (13) (if we perform the replacement z = in 
equation (13)). 
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3 Numerical Results 



The solutions (or the functions A{z) and 4>{z))) are obtained by solving numerically the 
system of second order differential equations (15) and (16). The boundary conditions on the 
brane 0(0), A'(0 + ) and 4>'(0 + ) are obtained by solving numerically the algebraic system of 
Eqs. (17), (18) and (13). Note that A(0) = 0, as a(0) is normalized to unity, and a(z) = e A{ - z \ 

We would like to emphasize that the parameter space of the model includes five param- 
eters £, A, A, a, b. Before we proceed with the numerical analysis, it is convenient to make 
the parameters dimensionless by performing the transformation z — > kz. In order to use 
the notation of the RS-model, without the scalar field, we choose k = ^J—A/6 (note that 
in this paper we have examined mainly the case of a negative five dimensional cosmological 
constant, or A < 0). Now the number of independent parameters is reduced to four, namely: 
£, A = A/A; 2 , a = cr/k, b = bk 2 . After this, it is obvious that the relevant scale is the bulk 
cosmological constant. 

Of course, a complete investigation of a so extensive parameter space, with four indepen- 
dent parameters, is beyond the scope of this paper. Thus, we will restrict our study to the 
classes of numerical solutions we obtained in j9]. We aim to investigate how Gauss-Bonnet 
gravity affects these solutions. 



3.1 Conformal Coupling 




Figure 1: The warp factor a(z) and the scalar field 0( function of kz, for £ = £ c , 

a = 5, A = 0.01 for two values of b: b = —0.1 (upper panel) and b = 0.04 (down panel). The 
values for a and A do not satisfy the fine tuning of the RS-model. 

In the case of 5D conformal coupling £ c = 3/16 and b = 0, as we have already mentioned 
in [9], the warp factor a(z) is of the order of unity in a small region near the brane (starts 
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with a'(0) negative) and increases exponentially (a(z) — » e kz ), as z — > +00, or equivalent the 
space-time is AdS$ asymptotically, while the scalar field <fi(z) is nonzero on the brane (in the 
region where the warp factor is almost constant) and tends rapidly to zero for large z. 

In Fig. [1] we see that the above mentioned behavior for the warp factor and the scalar 
field remains, even when the Gauss-Bonnet term is present. In particular we have plotted 
the warp factor a(z) and the scalar field <p( function of z, for £ = £ c , for two values 

of b: b= -0.1 and b = 0.04. 

A complete investigation of the spectrum of the solutions of the model is quite com- 
plicated, however we will try to give a description. We have checked, that for b < 0, the 
algebraic system of Eqs. (17) and (18) (junction conditions) has only one real solution, which 
corresponds to a warp factor and a scalar field profile, of the form which is plotted in the 
upper panel of Fig. [TJ On the other hand, for < b < 6 c3 (6 c3 = 0.16), the algebraic system 
of Eqs. (17) and (18) has three real solutions which correspond to three static solutions of 
the Einstein equations. For < b < b c \ (b c i = 0.1227) the two of them are of the form 
of the down panel of Fig. [T] (Ads?, in the bulk), and the third has a warp factor which 
vanishes in finite proper distance z and <fi(z) tends to infinity there (or equivalently it has 
a naked singularity in the bulk). In the short interval b c \ < b < b C 2 (b C 2 = 0.1250) only one 
of the solutions is of the form of the down panel of Fig. [H and the other two exhibit naked 
singularities. For b C 2 < b < 0.16 we have three solutions with naked singularities. Finally, 
for b > 0.16 the algebraic system of Eqs. (17) and (18) has only one real solution, which 
exhibits a naked singularity in the bulk. 

The above presented analysis is valid for a special choice of the free parameters of the 
model (£ = £ c , a = 5, A = 0.01). For a different choice the analysis of the spectrum of 
the solutions is quite similar. For example, in the case of (£ = £ c , a = 2, A = 0.01), we 
obtain that b c \ = 0.1247 and b C 2 = 0.1250. Note that the second critical value b C 2 is almost 
independent from the value of brane tension a. For large b (b > b c ^, b c3 = 1 ) the model has 
only one solution with a naked singularity in the bulk. 

In the case of a zero tension brane a = 0, A = 0.01 we have also solutions of the form 
we described above, with concentrated around zero (z = 0) and the space-time is AdS§ 
asymptotically. For < b < 0.1250 we have three distinct solutions. Two of them are of the 
form of Fig. [T] and the third exhibits a naked singularity, with a zero warp factor in finite 
distance and infinite scalar field. Note that the first of the nonsingular solutions is smooth, 
while the second exhibits an effectively negative brane tension due to the discontinuous first 
derivative of the scalar field for z = 0. However, if b exceeds the critical value b c = 0.1250 
we obtain three solutions which suffer from naked singularities in finite proper distance in 
the bulk with infinite Ricci scalar. Note that the two critical values b c \ and b C 2 which are 
different in the case of the nonzero tension coincide in the case of a zero tension brane a = 0, 
also fe C 3 goes to infinity. 

3.2 Class (a) (f < 0) 

The solutions for £ < 0, as we have seen in our previous paper of Ref. [9], exhibit a naked 
singularity (in particular the warp factor becomes zero in finite proper distance in the bulk), 
while the scalar field tends to infinity near the singularity. 

It is believed that higher order gravity may incorporate quantum gravity effects which 
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Figure 2: The warp factor a(z), as a function of kz, for b = —0.1 , a = 6, A = 0.01, for two 
values of £: £ = —0.001 (left-hand panel) and £ = —0.1 (right-hand panel). The values for 
a and A satisfy the fine tuning of the RS-model. 



are important when a singularity appears in a solution of Einstein equations. The question 
that arises is whether Gauss-Bonnet gravity can resolve the naked singularities we obtained 
in p. 

In order to answer the above question, we examined a wide range of the parameter space 
of the model, with £ < 0, and we found that Gauss-Bonnet gravity can not overcome the 
singularity in a satisfactory way. 

A characteristic behavior for the solutions is presented in Fig. [5J where we have plotted 
the warp factor a(z), as a function of kz, for b = —0.1 , a — 6, A = 0.01, and for two 
values of £: £ = —0.001 (left-hand panel) and £ = —0.1 (right-hand panel). We see that 
naked singularity in the bulk remains, even when the coupling b possesses a nonzero value. 
Note that the scalar field tends to infinity near the naked singularity. The warp factor for 
b = is presented in Figs. (1) and (2) in Ref. [9]. If we compare with the case of nonzero 
b of Fig. [2] in this work, we see that the singularity points are only slightly displaced. We 
have examined also other values of negative b, and we observed that for any negative b the 
solutions are of the form we described above. 

For positive b we see that the algebraic system of equations (13), (17) and (18) has one 
or more real solutions. Even in this case we cannot avoid naked singularities in the bulk. 
However, these new singularities are of a different nature from the one we described above. 
In particular in this case, the second derivative of the warp factor (or the Ricci scalar) tends 
to infinity in finite proper distance in the bulk, while the warp factor, the scalar field and 
the first derivative of the warp factor tend to finite values. 

For zero tension branes a = we observe the same behavior with that described above. 
For b < we have one singular solution with a warp factor which tends to zero, while for 



7 



b > we have only one singular solution with nonzero warp factor but infinite Ricci scalar. 
3.3 Class (b) (f > f c ) 
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kz kz 



Figure 3: The warp factor a(z), as a function of kz, for £ = 1 + £ c , er = 6, A = 0.01, for 
two values of b: b = —0.1 (upper panel) and b = 0.01 (down panel). The values for a and A 
satisfy the fine tuning of the RS-model. 

For 6 = and £ > £ c , we had obtained a different class of numerical solutions in Ref. [9]. 
In this case the warp factor a(z) is of the order of unity in a small region near the brane 
and increases exponentially (a(z) — > e fez ), as 2 — > +oo (or equivalently the space-time is 
asymptotically AdS^), while the scalar field <J)(z) is nonzero on the brane (in a rather small 
region near the axes origin where the warp factor remains constant), and for kz » 1 tends 
rapidly to zero. 

The question that arises is, how these solutions are modified when a Gauss-Bonnet term 
is turned on. As we see in the upper panel of Fig. [3X for b = —0.1, the features of the solution 
remain the same with those for 6 = 0. In general, for b < 0, we have only one static solution 
of the form we described. We have checked numerically that for < b < b c \ (b c i = 0.0106) 
the model possesses two static solutions which are asymptotically AdS$. The first of them 
is of the previously described form, while the second has the same profile with the first, but 
it is translated to a finite proper distance in the bulk, at it is presented in the down panel 
of Fig. [21 For b cl < b < b c2 (b c2 = 0.0192) we have also two solutions, but only one is of the 
form we described above, while the other exhibits a naked singularity. For b > 6 c2 we obtain 
one or more singular solutions depending on the value of b. 
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Figure 4: The warp factor a(z) and the scalar field <f>( function of kz, for £ = 0.1, 

a = 6, A = 0.01, for two values of b: b = 0.5 (upper panel) and b = —0.04 (down panel). In 
the case of b = 0.5 we have naked singularity at z s = 0.42 in the bulk. 

3.4 Class (c) (0 < £ < f c ) 

In the absence of the Gauss-Bonnet term (b = 0) and for < £ < £ c the warp factor a(z) 
and the scalar field <fi(z) tend rapidly to infinity near a singular point in the bulk, like the 
upper panel of Fig. 0] (for details see Ref. [9]). 

We will present an investigation of the spectrum of the static solutions of the model in 
the case of a specific choice of the free parameters of the model £ = 0.1, a = 6 and A = 0.01. 
For < b < b c i (b c \ = 0.0938) we have two solutions with singular behavior in the bulk. The 
first of them exhibits the behavior we described for 6 = 0, like the upper panel of Fig. HJ and 
the other has an infinite second derivative of the warp factor, while the functions a(z), a'(z) 
and (p(z) are finite near the singular point. For b > b c \ we obtain only one solution of the 
form the upper panel of Fig. HI In the case of negative b and b c2 < b < (b c2 = —0.0485) 
we obtain three distinct solutions: a) Singular solutions where the warp factor a(z) vanishes 
in finite proper distance in the bulk, b) Singular solutions where the second derivative of 
the warp factor tends to infinity, while the functions a(z), a'(z) and <j>{z) are finite near 
the singular point in the bulk and c) Nonsingular solutions with a warp factor a(z) which 
increases exponentially (a(z) — > e kz ), as z — > +oo and a scalar field <f)(z) which is nonzero 
on the brane, and for kz » 1 tends rapidly to zero, see the down panel of Fig. HI For 
b c s < b < b c2 (6 C 3 = —0.0536) we have three singular solutions, and for b < b c3 we have only 
one singular solution. 

In general, we observe a similar investigation for other values of the free parameters of 
the model. In particular, in the case of a = we obtain one singular solution of the form of 
the upper panel of Fig. HJ for all the values of the Gauss-Bonnet coupling b. 
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4 Conclusions and discussion 



We considered brane world models with an additional Gauss-Bonnet term in the presence 
of a massless bulk scalar field which interacts non-minimally with gravity, via a possible 
interaction term of the form — ^Rcj) 2 . We developed a numerical approach for the solution 
of the Einstein equations with the scalar field, and we studied the spectrum of the static 
solutions, in the case of a A0 4 potential, for several values of the free parameters of the model. 
We obtained that the existence of the Gauss-Bonnet term does not change significantly the 
solutions we found in our previous work [5J, in appropriate regions of the Gauss-Bonnet 
coupling b. Outside these regions the solutions always suffers from naked singularities at 
finite proper distance in the bulk. In general the spectrum of the solutions is characterized 
from naked singularities except for three cases: i) £ > £ c and b < b c (unstable solutions 
against scalar field perturbations with a tachyonic spectrum), ii) £ = £ c and b < b c (stable 
solutions against scalar field perturbations with a continuous spectrum of positive energies) 
and iii) < £ < £ c and b c < b < (stable solutions against scalar field perturbations 
with a discrete spectrum of positive energies). Note, that the third case of stable solutions 
is new and does not appear when b is equal to zero. Furthermore these solutions share 
the same properties: a(z) increases exponentially with z, and tends to an AdS$ space- 
time asymptotically, while (f)(z) is nonzero near the brane and vanishes rapidly when z is 
increasing. The brane at z = has a positive (or zero) brane tension, and the cosmological 
constant in the bulk is negative. We note, that in the presence of the Gauss-Bonnet term 
we obtain solutions which are AdS$ asymptotically even with positive bulk cosmological 
constant, but we have not present our study in this case. 

Finally, it would be interesting to discuss if the solutions we obtained can be used for the 
construction of a more realistic brane world model. The main motivation to consider these 
models was the observation that the interaction term — |£i?0 2 can be interpreted as a mass 
term, with an effective mass m 2 ff = £R. In the case of the second RS-model we obtain that 
R = 16k5(z) — 20k 2 , and we see, that for negative £, the effective mass term is negative on 
the brane and positive in the bulk. This implies an instability of the RS-vacuum near the 
brane. As we argued in Ref. [10J, we expected static stable solutions with a warp factor 
similar with that of RS-model, and a scalar field vacuum with nonzero value on the brane 
and zero value in the bulk. These kind of solutions incorporate a layered-phase mechanism 
for the localization of Standard-Model particles on the brane. This mechanism is based on 
the idea of the construction of a gauge field model which exhibits a non-confinement phase 
on the brane and a confinement phase on the bulk. Thus gauge fields, and more generally 
fermions and bosons with gauge charge, can not escape into the bulk unless we give them 
energy greater than the mass gap Ac, which emerges from the nonperturbative confining 
dynamics of the gauge field model in the bulk. For an extensive discussion see |10j . 

In [9J we solved the Einstein equations and we found that there not solutions with the 
profile we described above for £ < 0. However we obtained a class of stable solution, for 
£ — £c which may be physically interesting. In particular, the warp factor a(z), for this class 
of solutions, is of the order of unity near the brane and increases exponentially (a(z) ~ e fc ' 2 
), as z — > ±oo, while the scalar field <j>{z) is nonzero on the brane and tends rapidly to 
zero in the bulk. The model is completed if we include a second positive tension brane with 
a'/k = 6, in a position z c where the scalar field vacuum is practically zero. This two-brane 



10 



set up is very similar to the well known RSI-model. In analogy with the RSI-model, we will 
assume that the standard model particles live in the effectively negative tension brane 0, or 
visible brane. The advantage of this model is that it incorporates a layered-phase mechanism 
for the localization of standard model particles on the visible brane, and analyzed in detail 
inRef. 0. 

In this paper we generalize our previous works by considering brane-models with an 
additional Gauss-Bonnet term and a bulk scalar field nonminimally coupled with gravity. 
As we have already mentioned the addition of the Gauss-Bonnet term has not a significant 
impact to the spectrum of the model, and gives solutions with the same features as the model 
without the Gauss-Bonnet term. However this model exhibits solutions of the interesting 
form we described in the previous paragraph, even when £ 7^ £ c . In particular these solutions 
appear when < £ < £ c and b c < b < 0, and can be used for the construction of realistic 
brane world models, as we describe above. 

We would like to note that in the phenomenological studies in [16] the same model, 
with an additional Gauss-Bonnet term and a bulk scalar field nonminimally coupled with 
gravity, is considered for £ < 0, however the authors of this paper do not take into account 
the back reaction of the scalar field to the RS vacuum. We think that is important for 
the phenomenological applications to have a stable background solution, for the complete 
Einstein equations and the scalar field. 
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